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We present a search for the radiative leptonic decay B + — » £ + vi~/, where I = e, /J,, using a data 
sample of 465 x 10 6 BB pairs collected by the BABAR experiment. In this analysis, we fully re- 
construct the hadronic decay of one of the B mesons in Y(4S) — ► B + B~ decays, then search for 
evidence of B + —> l^va in the rest of the event. We observe no significant evidence of signal decays 
and report model-independent branching fraction upper limits of B(B + — > e + v e y) < 17 x 1CF 6 , 
B{B + ~* < 24 x 1(T 6 , and B(B + -> < 15.6 x 1(T 6 [I = e or n), all at the 90% 

confidence level. 
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The leptonic decay i + v^ where I = e or p,, 

proceeds via quark annihilation into a virtual W + bo- 
son with the radiation of a photon. The presence of the 
photon removes the helicity suppression of the purely 
leptonic decays, B + — > £ + vi, although it introduces an 
additional suppression by a factor of a cm . The branch- 
ing fraction of B + — ► l + v^ is predicted to be of order 
10 -6 [2j, making it potentially accessible at B factories. 
The most stringent published limits are from the CLEO 
Collaboration with B{B+ -> e+v e j) < 2.0 x 1(T 4 and 
B(B+ -> n + Vu-f) < 5.2 x 10~ 5 at the 90% confidence 
level (C.L.) |jf. 

The differential branching fraction versus photon en- 
ergy E~f involves two form factors, fv and /a, which 
contain the long-distance contribution of the vector and 
axial currents, respectively, in the B — > 7 transition 



dB _ a enl G F \V ub \ 2 



48tt 2 



m%r B [f%{E^) + f v (E 7 j] (l-y)y 3 , 

(1) 

where Gf is the Fermi constant, V u b is the Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix element de- 
scribing the coupling of b and u quarks, raj and tb 
are the JS-meson mass and lifetime, respectively, and 
y = 2E~ l /mB- While Ja — fv m most models [4], some 
suggest $a = [1| . The branching fraction is given by 
Ref. S as 



B{B+^ l+vef) 



Qcm<j|.|V^i,| 
2887T 2 



., 1 Qu Qb 

.f B m B T B t 



A, 



(2) 

where fg is the _B-meson decay constant, Q u ,b are the u- 
and 6-quark charges, and rat is the 6-quark mass. The 
first inverse moment of the B-meson distribution ampli- 
tude As is expected to be of order Aqqu but its theoreti- 
cal estimation suffers from large uncertainties 0] • It also 
appears in the branching fractions of two-body hadronic 
-B-meson decays, such as B — > it-it, and plays an impor- 
tant role in QCD factorization [4j. Since there are no 
hadrons in the final state, an experimental measurement 
of B + — > £ + vi^/ can provide a clean determination of As. 

We present the first search for B + — > t^vti that ex- 
ploits the hadronic "recoil" technique, in which one B 
meson is exclusively reconstructed in a hadronic final 
state before searching for the signal decay within the 
rest of the event. This technique improves the handling 
of event kinematics, providing adequate background sup- 
pression without requiring model-dependent constraints 
on the signal kinematics. Thus, this analysis is valid for 
all B — > 7 form-factor models and over the full kine- 
matic range. This analysis uses a data sample of 465 ± 5 
million BB pairs, corresponding to an integrated lu- 
minosity of 423 ftT 1 collected at the T(45) resonance. 
The data were recorded with the BABAR detector at the 
asymmetric-energy PEP-II e + e~ storage ring at SLAC. 
The BABAR detector is described in detail elsewhere [8J . 



Signal and background decays are studied using Monte 
Carlo (MC) samples based on GEANT4 [§]. The sim- 
ulation includes a detailed model of the BABAR detec- 
tor geometry and response. Beam-related background 
and detector noise are extracted from data and over- 
laid on the MC simulations. 7(45) -> BB signal MC 
samples are generated with one B meson decaying via 
B + ^ l^vfy using the tree-level model of Ref. [6j, which 
is valid for y > 0.13, while the other B meson decays 
generically. We simulate signal MC samples for two form- 
factor models, with f a = fv arid Ja = respectively, to 
evaluate the impact of the decay model on the signal 
selection efficiency. Large MC samples of generic BB 
and continuum (e + e~ — > t + t~ or e + e~ — > qq, where 
q = u,d,s,c) events are used to optimize the signal 
selection criteria. However, the final background esti- 
mates are obtained directly from a combination of data 
and exclusive B + — > X^£ + vi MC samples, where X° is 
a neutral meson containing a u quark. The primary 
background for B + — > in this analysis is due to 

B+^ X°£+v t decays, with B+^> ir°l + v t {B+^ r)£+v e ) 
comprising approximately 73% (18%) of this semileptonic 
background. The branching fraction and uncertainty for 
each B + — > X^t + V£ mode are taken from experimental 
measurements (X° = it 0, p 0, v EH and to 0). 
We assume B(B+^ r)'£+v e ) = B(B+^ rji+vg) x (1 ± 1). 
We use a light-cone sum rule model for the 77 and rf form 
factors [13| and use the form factor measured in a BABAR 
analysis [14| with the shape parameterization given in 
Ref. [HI for the tt° mode. 

Event selection begins with the full reconstruction of 
a charged B meson (-Btag) in one of the large number 
of hadronic final states, B~ — » -D^Xhad- We recon- 
struct the D*~ — > D°tt-; D*° -> D tt°, D°r, D~ -> 

A'°7T~, K®1T~1T , K° S 1T~ 1T~ 1T + , K + 7T~1T~, K + 7T~ 7T~ TT° ] 

D° -> K+tt-, K+ir-iT , K + it-it-it+, K°tt-tt+; and 
K®-^ it~tt + decay modes. Xhad is a collection of at most 
five mesons, composed of both charged and neutral kaons 
and pions. Well- reconstructed -B tag candidates are se- 
lected using two kinematic variables: AE — Es t&g — \fs /2 

and m ES = - p Bt ^, where E Btas and p Btas are 

the energy and momentum of the Stag candidate, re- 
spectively, and \/s is the total energy of the e + e~ sys- 
tem, all in the center-of-mass (CM) frame. We require 
AE, which peaks at zero for correctly reconstructed B 
mesons, to lie between —0.12 and 0.12 GeV or within two 
standard-deviations from its mean for the given Ah a d 
mode, whichever is the tighter constraint. We fit the 
toes distribution for each Ah a d mode and require that 
the purity, or fraction of well-reconstructed B mesons, is 
greater than 12% in the region toes > 5.27GeV/c 2 . If 
more than one -Btag candidate is reconstructed, the one 
in the highest purity mode is chosen. If there are multi- 
ple candidates in this mode, the one that minimizes \AE\ 
is selected. 
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FIG. 1: mES distribution, after Stag reconstruction and con- 
tinuum suppression, of data (points) and the expected com- 
binatoric background as predicted by the MC (shaded). 



We define the signal region as 5.27 < m ES < 5.29 
GeV/c 2 , since correctly reconstructed B mesons peak in 
this region near the nominal i3-meson mass. The B tag 
candidates that are incorrectly reconstructed from ei- 
ther continuum events or both B mesons ( "combinatoric" 
events), produce a distribution that is fairly flat below 
the ttt-es signal region and decreases within it, as shown 
in Fig. [T] The shape of the combinatoric distribution is 
extrapolated into the toes signal region using MC, while 
the background contribution from combinatoric events is 
estimated directly from the data. To improve the MC es- 
timate of the Stag reconstruction efficiency, we normalize 
the generic MC to the number of data events that peak 
within the toes signal region. Thus, all MC samples are 
scaled by 90.7%, resulting in good agreement between 
data and background MC throughout the analysis selec- 
tion. A charged B tag is reconstructed in about 0.3% of 
the signal MC events. 

Because the two B mesons produced in the T(4S) de- 
cay have low momenta in the CM frame (0.3 GeV/c), their 
decay products are more isotropic than continuum back- 
ground. For example, |cos#t|, where 6t is the angle 
in the CM frame between the B tag thrust axis and the 
thrust axis of all other particles in the event, has a flat 
distribution for BB events and peaks near one for non- 
BB events. The continuum background is suppressed by 

requiring C B = 11^0*)/ (IL^uO*) +IL7V*i)) > 
30%, where Vsixi) (V q {xi)) are probability density func- 
tions determined from MC that describe BB (contin- 
uum) events for the five event-shape variables x\. The 
variables used are: the ratio of the second to zeroth Fox- 
Wolfram moment [Tsl ] computed using all charged and 
neutral particles in the event, the cosine of the angle be- 
tween pB tas and the beam axis, the magnitude of the J5 tag 
thrust, the component of the Stag thrust along the beam 
axis, and |cos0t|- This requirement improves the agree- 
ment between data and MC by suppressing unmodeled 
continuum backgrounds, such as e + e~— > e + e~£ + £~ via 
two photons. 
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FIG. 2: m„ distribution after all selection criteria are ap- 
plied, in electron (top) and muon (bottom) modes for the 
WES-peaking (shaded) plus nonpeaking (solid) contributions 
in the full background MC sample, signal MC normalized to 
B — 40 x 10~ 6 (dashed), and data (points). Events to the 
left of the vertical lines are selected. iV f comb of Table U is de- 
termined from sideband data, not from the MC shown here. 



In the sample of selected B tag candidates, we iden- 
tify events in which the remaining tracks, calorimeter 
clusters, and missing momentum vector (p m i ss ) are con- 
sistent with B + — * l + vi^ candidates. We select events 
with exactly one track, which reduces the signal effi- 
ciency by 25% but removes over 99% of the simulated 
background events with a reconstructed -Btag- This sig- 
nal track is required to have a charge opposite to that 
of the Stag, to satisfy particle identification (PID) crite- 
ria for either a muon or electron, and to be inconsistent 
with a kaon hypothesis. In the electron mode, the four- 
momenta of signal tracks are redefined to include those of 
any bremsstrahlung photon candidates. Such a candidate 
is defined as any cluster whose momentum vector, when 
compared to that of the signal track (pi), is separated by 
|A0| < 3° and -3° < Q e x Acp < 13°, where Q e = ±1 
is the e charge and 9 (4>) is the polar (azimuthal) an- 
gle relative to the beam axis, in the lab frame. Finally, 
the signal photon candidate is chosen as the cluster with 
the highest CM energy, excepting bremsstrahlung photon 
candidates. 

We significantly reduce the background by requiring 
that the kinematics of the signal track and photon can- 
didate are consistent with the existence of a third mass- 
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less particle originating from the signal B meson. To do 
this, we use the four-momentum of the expected signal 
B meson (pb), which is assigned an energy of y/s/2, a 
momentum vector pointing along — Ps tag , and the nom- 
inal B-meson mass. The neutrino mass squared is then 
defined as m 2 = (pB—pi—p-y) 2 , where pi (j> 7 ) is the four- 
momentum of the signal track (photon candidate). As 
shown in Fig. [2j the background increases with m 2 , while 
B + ^r £ + i/^ events peak at m 2 = with an enhanced tail 
in the electron mode due to unrccovcrcd bremsstruhlung 
photons. We require -1 < m 2 < 0.46 (0.41) GeV 2 /c 4 
for the electron (muon) modes. In addition, the lepton 
and neutrino should be emitted back-to-back in the rest 
frame that recoils from the photon emission, defined as 
PB~P-y- We require cos0g„ < —0.93 in this frame, where 
6(, v is the angle between pi and p m i SS . After all other 
selection criteria are applied, the MC indicates that m 2 
and cos 9^ together remove 99% of background events 
with a 30 and 20% reduction in the signal efficiency for 
the electron and muon modes, respectively. 

The dominant backgrounds are due to B + ^ ir°£ + i/£ 
(r]£ + vi) events in which 7r°(?7)^ 77 fakes the B + ^ £ + vij 
signal photon. To suppress this background, we reject 
events containing a ir°(r]) candidate, reconstructed using 
the signal photon candidate and a second cluster hav- 
ing CM energy E l2 . For 7r° candidates, we require a 
77 invariant mass between 120-145 MeV/c 2 with E l2 > 
30 MeV or between 100-160 MeV/c 2 with E l2 > 80 MeV. 
For rj candidates, we require a 77 invariant mass be- 
tween 515-570 MeV/c 2 with E l2 > 100 MeV. Likewise, 
B + ^ uj£ + i>t—> [7r°7]£ + ^ events are suppressed by reject- 
ing any event in which the signal photon candidate and a 
7T° candidate produce an invariant mass between 730-830 
MeV/ c 2 . This ir° candidate is defined as any two clusters 
with CM energy > 70 MeV which produce a 77 invariant 
mass between 115-145 MeV/c 2 . After applying all other 
selection criteria, these vetoes reduce the B + ^> n°£ + ve. 
and B + — > X®£ + vi background events, with AT" 7^ 7r°, 
by 65% and 50% respectively. Finally, we require the 
lateral moment [13] of the calorimeter energy deposit for 
the signal photon candidate, which peaks at 25% for sin- 
gle photons, to be between and 55%. This suppresses 
B + ^ ir°£ + i>i events in which the two photons from the 
7T° decay are reconstructed as a single merged photon. 

Once the B tag , signal photon, and lepton are identified, 
B + ^r £ + vgj events are expected to contain little or no 
additional energy within the calorimeter. However, ad- 
ditional energy deposits can result from hadronic shower 
fragments, beam-related photons, and photons from un- 
reconstructed D*— ► Dj/tt° transitions in the Stag can- 
didate. The total energy of all additional clusters is re- 
quired to be less than 0.8 GeV, counting only clusters 
with lab-frame energy greater than 50 MeV. We also re- 
quire that Pmiss points within the fiducial acceptance of 
the detector. 

To avoid experimenter bias, we optimize all the se- 




m 2 v for B + -> (GeV 2 /c 4 ) 

FIG. 3: distribution for B + ^> -k°£ + vi (I = e and fj,), using 
the procedure described in the text where 7 is substituted 
with a 7T° candidate, of data (points) and of B + — > ty°£ + ui 
MC normalized to B = 7.7 x 1CT 5 (dashed) and added to the 
expected background (solid). 

lection criteria and determine the number of expected 
background events in the signal region (A^ kg ), for £ = e 
or /i, before looking at any data events selected by the 
criteria. We optimize by maximizing the figure of merit 

4 g l(\ n ° + \f N e kS ) 03. where n ° = L3 and £ t S is thc 
total signal efficiency including that of the -Btag recon- 
struction. The signal branching fraction is calculated us- 
ing & = (N° bs -N bkg )/ef 8 N B ±, where N B ± = 465x 10 6 
is the number of B^ mesons in the data sample and A^ bs 
is the number of data events within the signal region. 

To verify the modeling of e^ Ig , we remove the B + — > 
X^£ + vi vetoes, select events containing a ir° candidate, 
and substitute the 7r° in place of the signal photon can- 
didate. The resulting m 2 distribution from B + ^ ir°£ + V( 
is expected to resemble that of the signal. We observe 
a peak in the data that agrees with MC expectations 
within the 15% statistical uncertainty of the data, as 
shown in Fig. [3l For cross-check purposes only, we de- 
termine the B + — ► ir°£ + i>i efficiency using an exclusive 
B + ^> n°£ + vi MC sample and the background contribu- 
tion using generic MC. The peak in data corresponds to 
B(B+^ n°£ + ve) = (7.8±^) x 10" 5 , where the uncer- 
tainty is statistical. This branching fraction is consistent 
with the current world-average value of (7.7± 1.2) x 10~ 5 
[Io| . which is also the value used in the MC samples. 

The total number of background events N^ ks has two 
components: Nf e the number of expected background 
events having a correctly reconstructed Stag and hence 
peaking within the mgs signal region, and A^ omb the 
number of expected combinatoric background events, in- 
cluding both BB and continuum events. The m 2 and 
cos 6i v restrictions ensure kinematic and topological con- 
sistency with a three-body decay involving a massless and 
undetected particle: the neutrino. By further requiring 
that exactly one track recoils from a fully-reconstructed 
^tag, lepton number and PID ensures the track is a lep- 
ton. Thus, only B + — > £ + V(j decays can peak within 
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the signal region, unless the signal photon candidate ac- 
tually arises from one or more particles that mimic the 
kinematics of B + — > £ + is£j, which only occurs in spe- 
cific pathological B + — > X®£ + i/£ decays. Therefore, we 
determine Nf eak using exclusive X®£ + vi MC sim- 

ulations and validate the lack of additional peaking back- 
grounds with generic MC. Other decay modes passing the 
selection criteria do so with poorly reconstructed Stag 
candidates and thus produce a combinatoric distribution 
in TTiEs- We determine iV™ mb from an extrapolation of 
the observed number of data events within the toes side- 
band region, defined as 5.20 < Toes < 5.26GeV/c 2 . We 
observe 1 (4) data events within the Toes sideband for 
the electron (muon) mode. 

The uncertainty on iV™ mb is dominated by the side- 
band data statistics. It also includes the systematic un- 
certainty from the combinatoric background shape, esti- 
mated by varying the selection criteria and the method 
used to extrapolate this shape (14.6%). The error on 
Nf eak is dominated by uncertainties in the branching 
fractions and form factors associated with the various 
exclusive X®£ + vi decays (13.6%). Additional sys- 

tematic uncertainties result from MC modeling of the 
data efficiency, which we apply to both Nf ea and e^ lg : 
electron PID (0.9%) or muon PID (1.3%), C B (1.4%), to 2 
(0.5% for e s e ls , 1.4% for Nf cak ), and the reconstructions 
of the track (0.4%), photon (1.8%), and B tag (3.1%). The 
last of these, which also accounts for uncertainty in N B ± , 
is estimated by varying the shape of the toes combina- 
toric distribution and the size of the toes signal and side- 
band regions. 

Branching fraction limits and uncertainties are com- 
puted using the frequentist formalism of Feldman and 
Cousins [13], with the uncertainties on N^ ks and e^ Ig 
modeled using Gaussian distributions. Since B(B + — > 
£ + vi"f) is expected to be independent of the lepton type, 
we also combine the two modes by maximizing a like- 
lihood function defined as the product of both Poisson 
probabilities in N^ kg , where Bg is the mean. 

We observe 4 (7) data events within the signal region 
for the electron (muon) mode, compared to an expected 
background of 2.7 ± 0.6 (3.4 ± 0.9) events. This cor- 
responds to a signal significance of 1.2er (1.8c), a com- 
bined significance of 2.1er, and the results given in Ta- 
ble Q] The effective detector and PID thresholds are 
about 20MeV for photon energy and 400 (800) MeV/c 
for electron (muon) momentum, and we apply no min- 
imum energy requirements. Thus, this analysis is es- 
sentially independent of the kinematic model; we as- 
sume the Ja = fv signal model, but the f A —0 model 
yields consistent e^ lg values. Since certain theoretical 
calculations are most reliable at high E 1 Q, we also 
report a partial branching fraction limit AB by select- 
ing events with a photon candidate energy greater than 
1 GeV, which reduces e s e ls by 30%. We observe 2 (4) data 



events with iV £ g = 1.4 ± 0.3 (2.5 ± 1.0), resulting in 
AB(B+^ l+vti) < 14 x 10- 6 at 90% C.L. 

In Table[H we also report model-specific limits by intro- 
ducing a kinematic requirement on the relationship be- 
tween cos# 7 £ and cos# 7 „, where 9-yg (6y U ) is the angle be- 
tween the photon candidate momentum and pg (pmiss) in 
the signal B rest frame. The photon is emitted preferen- 
tially back-to-back with the lepton in the Ja — Iv model, 
and with either the lepton or neutrino when Ja = 0. 
Thus, we require (cos9 7 e - l) 2 + (cos# 7 „ + l) 2 /3 > 0.4 or 
(cos6» 7l/ -l) 2 + (cos6» 7<! + l) 2 /3 > 0.4for the j a — model, 
and only the former relationship for f a = fv ■ This re- 
duces e s e ls in both modes and models by 40%. We observe 
(0) data events in the electron (muon) mode with A^" kg 
= 0.6 ± 0.1 (1.0 ± 0.4) for the f A = fv model, and 3 (2) 
data events with iV £ bkg = 1.2 ± 0.4 (1.5 ± 0.6) for f A = 0. 

In conclusion, we have searched for B + — > £ + vg-f us- 
ing a hadronic recoil technique and observe no signifi- 
cant signal within a data sample of 465 x 10 6 BB pairs. 
We present model-specific branching fraction limits in 
Table [T] We also report a model-independent limit of 
B(B+-+ £ + va) < 15.6 x 10~ 6 at the 90% C.L., which 
is consistent with the standard model prediction and is 
the most stringent published upper limit to date. Us- 
ing Eq.© with f B = 0.216 ± 0.022GeV [U, m B = 
5.279 GeV/c 2 , t b = 1.638 ps,ro b = 4.20GeV/c 2 , and 
\V ub \ = (3.93 ±0.36) x 10~ 3 [10], the combined branch- 
ing fraction likelihood function corresponds to a limit of 
A B > 0.3 GeV at the 90% C.L. 
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TABLE I: Expected background yields N e g = JVf om + iVj ea , signal efficiencies e B e ls , number of observed data events Nf 



resulting branching fraction limits at 90% C.L., and the combined central value ^combined- 
presented. Uncertainties are given as statistical ± systematic. 



Model-specific limits are also 



jycomb 
jypeak 

,sig 

£ e 

Nf° s 

^combined 

Model-independent limits 
/a = fv limits 
Ja — limits 



0.3 ± 0.3 ± 0.1 
2.4 ± 0.3 ± 0.4 
2.7 ± 0.4 ± 0.4 
(7.8 ± 0.1 ± 0.3) xl0~ 4 
4 

< 17xl0~ 6 

< 8.4 x 10~ 6 

< 29 x 10~ 6 



1.2 ± 0.6 ± 0.6 
2.1 ± 0.3 ± 0.3 
3.4 ± 0.7 ± 0.7 
3.1 ± 0.1 ± 0.3)xl0~ 4 
7 

< 26xl0~ 6 

< 6.7 x 10~ 6 

< 22 x 10~ 6 



(6.5±I:?±S:|) x io- 

< 15.6 x 10~ 6 

< 3.0 x 10~ 6 

< 18 x 10~ 6 
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